Web-assisted tunneling in the kicked harmonic oscillator 
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We show that heating of harmonically trapped ions by periodic delta kicks is dramatically en- 
hanced at isolated values of the Lamb-Dicke parameter. At these values, quasienergy eigenstates 
localized on island structures undergo avoided crossings with extended web-states. 
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Controlling the state and the time evolution of quan- 
tum systems is one of the central themes of current re- 
search in experimental and theoretical atomic physics, 
quantum optics, and mesoscopics. Tailoring wave pack- 
ets in Rydberg systems jlj, producing single photons on 
demand 0, creating coherent superpositions of macro- 
scopic persistent-current states 3], and, most recently, 
the controlled production of multiparticle entanglement 
are prominent examples of what is often coined as 
"quantum state engineering" , such as to stress the al- 
most perfect control we have achieved on matter on the 
microscopic level. Whilst many of these schemes still 
rely on an analysis of the quantum dynamics in terms 
of some unperturbed basis of the quantum system under 
control, it has become clear during the last decade that 
generic features of strongly coupled ("complex") quan- 
tum systems allow for novel and often extremely robust 
strategies of quantum control. In such systems, studied 
in much detail in the area of quantum chaos, peculiar 
eigenstates emerge which exhibit unexpected (and un- 
expectedly robust) localization properties and dynam- 
ics. Most prominent examples thereof are nondisper- 
sive wave packets in periodically driven quantum sys- 
tems ^ , quantum resonances , and quantum accelera- 
tor modes 0- A considerable part of these "strong cou- 
pling" quantum control schemes relies on some underly- 
ing classical dynamics, which in general is mixed regular- 
chaotic, and precisely the rich structure of a mixed classi- 
cal phase space is at their very origin d, E| ■ ^ large class 
of these systems follows the Kolmogorov-Arnold-Moser 
(KAM) scenario i-^- regular phase space structures 
are destroyed gradually as the coupling strength is in- 
creased. Yet, there is another kind of classically chaotic 
dynamics, which goes under the name "non-KAM" chaos 
[ll| , where the phase space flow is fundamentally altered 
at arbitrarily weak perturbations. Recently, the con- 
sequences of such non-KAM transitions have been ob- 
served in electron transport in superlattices where 
enhanced transport accross the lattice was observed, due 
to the sudden (non-KAM) appearence and disappearence 
of unbounded stochastic web-structures in classical phase 
space, at well-defined values of some control parameter. 
Hence, enhanced transport was enforced by (controlled) 
abrupt changes in the underlying classical phase space 
structure. 



Here we expose a different pathway for the controlled 
enhancement of transport in another non-KAM system, 
the kicked harmonic oscillator 0, Q, 0, 0|, at fixed 
phase space structure, by simply tuning the effective 
value of h. We will see that, for a suitable choice of 
the initial oscillator state, dramatic enhancement of the 
energy absorption by the oscillator from the kicking field 
can be achieved due to avoided crossings of localized (reg- 
ular) with extended (web-like) eigenstates of the kicked 
systems. We argue that our scenario can be easily ob- 
served in state-of-the-art experiments on harmonically 
trapped, kicked cold ions or atoms. 

The classical Hamiltonian of the kicked harmonic os- 
cillator 0,01 describes a harmonically trapped (in ID, 
with trap frequency v) particle of mass m, subject to a 
one dimensional, spatially periodic potential (wave vec- 
tor k — 27r/A, modulation depth A) which is periodically 
switched on and off at integer multiples of the kicking 
period r: 
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Further inspection of the Hamiltonian reveals that the 
classical phase space structure is completely determined 
by the parameters 



K : 
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and 



27r 



(2) 



which define the stochasticity parameter and the ratio 
between kicking and oscillator period, respectively. This 
is a consequence of the form-invariance of the equations 
of motion derived from under transformation to the 
scaled position and momentum coordinates v = kx and 
u = kp/m. 

A specific property of the kicked harmonic oscillator 
is that its phase space is unbounded, thus allowing for 
transport to infinity (which, in the trapped ion scenario, 
is tantamount to unbounded heating). Furthermore, it 
exhibits peculiar symmetry properties determined by the 
ratio q — {2tt/v)/t of oscillator and kicking period: For 
integer q, in addition to confined regular islands, it dis- 
plays a stochastic web (reaching out to infinity in phase 
space) with crystal [q £ qc = {3, 4, 6}) or quasi-crystal 
symmetry [q ^ qc), along which the system diffuses for 
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suitable initial conditions. This web is characterized not 
only by its symmetry but also by its thickness that broad- 
ens (shrinks) as the value of the stochasticity parameter 
K increases (decreases). 

With the harmonic oscillator annihilation operator 
a = [v + u)/'2ri and its hermitean conjugate at derived 
from the scaled center of mass coordinates of the trapped 
particle, the quantum version of ^ reads 

oo 

H = hiya^a + hK {cos [ri{a + a))] ]^5{t- nr), (3) 

n=0 

with K = K/2rf, and the Lamb-Dicke parameter rj = 
kyjh/2mv. Note that the latter measures the ratio of 
the width of the harmonic oscillator ground state in units 
of the wave length of the kicking potential, and that its 
squared value plays the role of an effective Planck con- 
stant - which can be tuned easily in ion trap experiments. 

Since (0 is periodic in time, the time evolution of any 
initial state ji/jo) is given by the action of the one-cycle 
Floquct propagator 

ij\ej) = e-'"«^"e-**™4')('i+a*)] = e"^^ |ej), (4) 

with the quasi-energies and \ej) the associated Flo- 
quet eigenstates. Due to the i5-term in (0), the Floquet 
propagator factorizes into a free evolution and a kicking 
part, and inspection of the former shows that the above 
classical condition for the emergence of a web structure 
in classical phase space is also quantum mechanically dis- 
tinguished, since integer q is equivalent to a (fractional) 
revival condition for the free evolution right upon the 
subsequent kick. One therefore expects quantum signa- 
tures of web-sustained transport in classical phase space 
- which we shall now establish by combining dynamical 
and spectral information. 

Let us start with Fig. ^ where we plot the mean en- 
ergy of the kicked (quantum) particle, initially prepared 
in the harmonic oscillator ground state, as a function of 
the number of kicks imparted on it by the external field, 
for slightly different values of the Lamb-Dicke parameter 
77, and for crystal symmetry (7 = 6. We see that a change 
of T] by merely one percent dramatically affects the en- 
ergy absorption by the particle from the field, eventually 
leading to quantum transport which is much more effi- 
cient than the classical excitation process. Note that the 
latter is unaffected by such tiny changes of 77, since the 
results displayed in the figure are obtained for fixed K 
and g, hence for fixed phase space structure! 

Whilst this latter observation appears puzzling on a 
first glance, it is readily resolved by inspection of the evo- 
lution of the quasienergy spectrum of Q as a function of 
•q. A global view thereof is provided by Fig. |21 for qua- 
sicrystal {q = 5) and crystal {q ~ 6) symmetries, respec- 
tively, and otherwise the same parameters as in Fig. ^ 
Only those eigenvalues associated with eigenstates which 




200 400 600 

Number of kicks 



FIG. 1: Mean energy as a function of the number of kicks, 
for q = G, K = 2.0, \'(po) — |0), and Lamb-Dicke parame- 
ter rj = 0.459, 0.464 and 0.469, at fixed classical phase space 
structure (eq. 0)- Also the classical time evolution is shown, 
for comparison. The inset zooms into a given region of the 
level dynamics in the lower panel of Fig. |21 (denoted by an 
arrow there). Circles and crosses indicate eigenstates with an 
overlap with the initial state \ipo) larger or smaller than 10~^, 
respectively. Dramatic enhancement of the energy absorp- 
tion by the trapped particle from the kicking field is observed 
exactly at the center (77 — 0.464) of the avoided crossing. 

have an overlap larger than 10^"^ with \ipo) are shown. 
For both symmetries, two distinct initial states {ipo) are 
considered: l-^o) = |0) (circles), and |0) displaced to the 
points (1.2,2.0) and (1.3,3.0), located in the chaotic re- 
gion of phase space, for q = 6 and q = b (dots), respec- 
tively. 

We see that, for the vacuum initial state, the level dy- 
namics are rather similar in the quasicrystal and in the 
crystal case, what is just the quantum signature of the 
locally similar phase space structures in the vicinity of 
the origin. In contrast, they are quite distinct for dis- 
placed initial conditions, since, under these premises, the 
quantum dynamics probes the rather distinct web struc- 
tures of the underlying classical phase space (illustrated 
by the insets in the figure, which show how a single tra- 
jectory explores phase space in the long time limit - note 
the different scales!). Independently of symmetry and 
initial condition, an abundance of avoided crossings is 
apparent from both plots - a typical signature of the 
strong coupling between the driving field and the center 
of mass motion of the trapped particle. Also note that, 
in the crystal case with in the chaotic domain, some 
avoided crossings are organized in a regular manner, what 
strongly suggests a semiclassical origin - which, however, 
is not relevant for our present purpose. 

For the vacuum initial state, both symmetries exhibit 
isolated or overlapping avoided crossings (of variable size) 
at well defined values of 77, which turn out to cause the 
enhanced quantum transport observed in Fig.^ The in- 
set in Fig. ^ zooms into the level structure around the 
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FIG. 2: Spectrum of the Floquet operator Q as a function 
of the Lamb-Dicke parameter rj, for K = 2.0, q = 5 (top) and 
q = 6 (bottom). Only eigenstates with an overlap larger than 
10~^ with the initial state lipo) are represented. Filled circles 
represent [i/jq) = |0), while dots refer to a displaced vacuum 
state centered at (1.3,3.0) (top) and (1.2,2.0) (bottom). The 
insets show the classical phase space explored by a single tra- 
jectory after 40000 kicks, when launched at a point near the 
origin, for the same K and q. The arrow in the bottom plot 
refers to the inset of Fig. 



avoided crossing at (r/ = 0.464; (f) — 1.35) indicated by an 
arrow in the lower panel of Fig. |2| Here, two quasiener- 
gies (as a matter of fact, two quasienergy bands [l^. 
due to the high spectral degeneracy induced by the sym- 
metries of the Hamiltonian) undergo an isolated avoided 
crossing, with maximal tunneling coupling between the 
associated eigenstates at its center at -q = 0.464. Note 
from Fig.^that enhanced quantum transport is observed 
precisely at this value. In contrast, for only slightly 
smaller or larger values of the Lamb-Dicke parameter, i.e. 
slightly shifted with respect to the center of the avoided 
crossing, the trapped particle absorbs energy with a much 
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FIG. 3: (Color online) Husimi representations for the eigen- 
states associated with the points A, B, C and D of the avoided 
crossing in the inset of Fig. Q Phase space is spanned by the 
coordinates v/2ri and u/2'q, respectivley. Clearly, the eigen- 
states localized in the vicinity of the hyperbolic point at the 
origin (B and C) dominate the dynamics in the asymptotic 
part {ri = 0.459 and i] = 0.469) of the anticrossing in the inset 
of Fig. Q whilst extended states localized on the stochastic 
web (A and D) exhibit equal weight at the center of the cross- 
ing parameter window for enhanced quantum transport. 

smaller rate, rather similar to the classical energy ab- 
sorption. Hence, the width of isolated avoided crossings 
induces the sharp response of the system to tiny changes 
in rj, opening a narrow parameter window in which the 
system can efficiently tunnel from its initial state {ipo) 
into another eigenstate which mediates rapid transport. 

But which eigenstates of the kicked harmonic oscil- 
lator promote enhanced diffusion? Let us inspect the 
phase space projection of the anticrossing states in the 
inset of Fig. n "on the left" (labels A and C) and "on the 
right" (labels B and D) of the crossing's center, which are 
represented by their Husimi functions in Fig. O As we 
approach the crossing from small values of ry, the lower 
branch of the anticrossing in the inset of Fig. ^ has the 
largest overlap with I'i/'o)- The corresponding eigenstate 
for 77 = 0.459 is plotted in the lower left panel of Fig.|31 
and is clearly localized on the hyperbolic fixed point at 
the origin, and on its first replica on the crystal lattice. 
It anticrosses with the state (A) shown in the upper left 
panel of the figure. Clearly, this latter state is extended 
(note that it has large densities far away from the ori- 
gin), localized on the stochastic web, reaching out to the 
effective boundary of the phase space which is defined by 
the finite size of the basis which we use for the diago- 
nalization of U. Since such web states - by their very 
nature - can never be converged numerically, we care- 
fully checked that the associated eigenvalues saturate at 
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FIG. 4: Mean energy (left y-axis) after 600 kicks, as a function 
of the Lamb-Dicke parameter tj, and for fixed classical phase 
space structure, with K — 2.0 and q = 6. The signal is 
strongly enhanced at those values of 77 which define avoided 
crossings between localized and web-like eigenstates. 



a finite basis size, and that their localization properties 
within a finite phase space domain remain structurally 
invariant when increasing the numerical basis. On the 
other side of the avoided crossing, at = 0.469, the anti- 
crossing states have exchanged their localization proper- 
ties, as illustrated in the right panels of Fig. O and as to 
be expected from the usual anticrossing scenario: Now 
the upper branch of the crossing provides the strongest 
support for j'^o), and the associated eigenstate (B) is once 
again localized on the hyperbolic fixed point at the origin, 
whereas the lower branches' eigenstate (D) is localized on 
the web. 

Consequently, the web eigenstates of the kicked har- 
monic oscillator lend support for enhanced quantum 
transport as observed in Fig. ^ at sharply defined val- 
ues of the effective Planck constant (parametrized by the 
Lamb-Dicke parameter), where near-resonant tunneling 
from the hyperbolic fixed point onto the web becomes 
possible. Since the Lamb-Dicke parameter can be tuned 
by variation of the trap frequency u or of the kicking 
lattice's wavelength 27r/fc, such web-enhanced tunneling 
from the harmonic oscillator ground state will have a 
marked signature on the experimentally observed heat- 
ing rate vs. f], as illustrated in Fig.^ Note that each of 
the peaks in the mean energy signal (after a fixed num- 
ber of kicks) plotted vs. r] in this plot corresponds to 
anticrossings (or, for the broader resonance structures, 
to overlapping avoided crossings) between localized and 
extended states alike the one identified in Figs. Q and O 
at fixed phase space structure. 

Tuning the Lamb Dicke parameter (at easily variable 
interaction time!) thus allows for sensitive probing of 
the local density of states by an experimentally robust 
observable, and, vice versa, for efficient control of the 
diffusion properties of the complex quantum system un- 



der study. Also note that Fig. 0] is reminiscent of con- 
ductance fluctuations which occur in solids meso- 
scopic systems and strongly driven atoms - 

which bear experimentally accessible information on the 
spectral structure that supports the detected probability 
current. 
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